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ABSTRACT 

The stellar disks of many spiral galaxies are twice as large as generally thought. We use archival 
data from the Galaxy Evolution Explorer mission (GALEX) to quantify the statistical properties of 
young stellar clusters in the outer, extended disks of a sample of eleven nearby galaxies. We find an 
excess of sources between 1.25 and 2 optical radii, R25, for five of the galaxies, which statistically 
implies that at least a quarter of such galaxies have this cluster population (90% confidence level), 
and no significant statistical excess in the sample as a whole beyond 2 optical radii, even though one 
galaxy (M 83) individually shows such an excess. Although the excess is typically most pronounced 
for blue (FUV — NUV < 1, NUV < 25) sources, there is also an excess of sources with redder 
colors. Although from galaxy to galaxy the number of sources varies significantly, on average, the 
galaxies with such sources have 75 ± 10 blue sources at radii between 1.25 and 2i?25- In addition, the 
radial distribution is consistent with the extended dust emission observed in the far IR and with the 
properties of Ha sources, assuming a constant cluster formation rate over the last few hundred Myrs. 
All of these results suggest that the phenomenon of low-level star formation well outside the apparent 
optical edges of disks (R ~ R25 ) is common and long-lasting. 

Subject headings: galaxies: evolution, fundamental parameters, spiral, structure 



1. INTRODUCTION 

Simple questions about galaxies can be difficult 
to answer. For example, the question "How 
large are disk galaxies?" has many answers. One 
could cite results regarding the edges of optical 
disks (jvan der Kruit fc Sear Id ll981lL of the star for- 
mation threshold (iKennicuttJ Il989f ). of H I disks 
Roberts fe Haynei 1 1 9941) . or of dark matter halos 



Zaritsky et alJ Il993f ). Even within each of those spe- 
cific aspects of spiral galaxies there are ambiguities. For 
example, recent results have shown both that the stel- 
lar d i sks are often not sharp ly truncated (jPohlen et alJ 
l2002HPohlen fe Truiillo1l2006f) . that star formation exists 
beyond t he radius of the sharp decline in star fo rma- 
tion rate (jFerguson et alj|1998tff hilker et al.ll2005f) . and 
that the main ingredient for star formation, mo lecular 
gas, exists at large radii (|Braine fe Herpinll200l . The 
discovery of radially extended star formation in at least 
some spiral galaxies has implications for a range of topics, 
including the growth and development of disks, star for- 
mation in low-metallcity, low-density environment, and 
galactic dynamics. So far, the presentat ion of extended 
star formation, eithe r as H II reg i ons ( Ferguso n et al 
1998) or UV knots dThilker et alj 1200a iBianchi et al 



200a iGil de Paz et al.l I2005D . has focused on the dis 



covery of the phenomenon, and as such, has been pre- 
sented in either individual galaxies or samples of two 
or three systems. A statistical study was presented by 
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iBoissier et al.l (|2006l ). but while it focuses on integrated 
luminosity to reach interesting results regarding the star 
formation threshold and initial mass function, it does 
not treat the star clusters separately. Is there any sys- 
tematic behavior in extended disk star cluster formation 
that we can identify? In particular, what is the per- 
vasiveness and radial extent of the phenomenon? We 
use this study to lay some groundwork for our kinematic 
and c hemical abundance studies o f outer-disk H II re- 
gions ([Christlein fc Zaritsky 2007a. b; Her bert-Fort et all 

We use archival data from the GALEX (Galaxy Evolu- 
tion Explorer) mission to study the phenomenon of UV- 
selected knots in the outer regions (R > R25) of galaxies. 
Because our goal is to determine statistical properties 
valid for disk galaxies as a class, we will stack the data 
for a suitable sample of galaxies and examine the angu- 
lar cross-correlation function between the UV knots and 
the parent galaxy. This study will then provide a mea- 
sure of the excess number of UV knots in the vicinity 
of nearby galaxies and therefore on the prevalence and 
radial range of the extended star formation. In §2 we 
describe the data and our treatment of them. In §3 we 
describe the cross-correlation results and in §4 we discuss 
some implications. 

2. SAMPLE SELECTION AND GALEX DATA 

We select our sample of face-on, well-resolved, disk 
galaxies from the extensi ve catalog of GA L EX g alaxy 
observations published by IGil de Paz et al.l (|2006f ) . To 
select face-on galaxies, we define the inclination using the 
ratio of major, a, to minor axis, b, appropriate for pro- 
jected circular disks and then select galaxies with an in- 
ferred inclination less than 45° . To select galaxies of large 
angular extent that still lie well within a single GALEX 
pointing, so that we obtain a reliable measurement of 
the background source density from the same image, we 
retain only those galaxies with 2 < a < 10 arcmin. 
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Fig. 1. — Source identification in a subfield around NGC 628. Greyscale represents the GALEX NUV image (the subsection is 20 arcmin 
across). Small black circles denote detected sources. The large white circle denotes the excluded region in our analysis around NGC 628. 



The semi- major axis i s taken from the compilation of 
iGil de Paz et al.l (|2006f ) and refers to R25, which we also 
interchangeably call the optical radius. Finally, to select 
disk galaxies, we require < T-type < 7, where T-Type 
characterizes morphology as descri bed in the Third Ref- 
erenc e Catalog of Bright Galaxies (|de Vaucouleurs et al.l 
119911 ) and our criteria span morphological types Sa to Sd. 
These selection criteria result in a sample of 20 galaxies. 
We then require that both high quality far and near UV 
data exist in the MAST archive at the Space Telescope 
Science Institute. Lastly, we reject NGC 1068 because 
of the highly non-uniform background in the image that 
has most likely variably affected source detection. We 
arrive at the final list of 11 galaxies presented in Table 1 

and refer to this as our base sample. 

We refer the read er to iMartin et al.l J2005) and 
iMorrissey et al.l (|2005h for a description of the mission 
and its on-orbit performance. The images and object 
catalogs used here are drawn from the database made 
available by the GALEX team via the MAST archive. 
The archival catalogs contain source positions and pho- 
tometry in both the FUV (1350 — 1750 A) and NUV 
(1750 — 27 50 A) channels and wer e constructed using 
SExtractor (|Bertin fc Arnoutdll996T ). Outside ~ 1.5 op- 
tical radii from the centers of the target galaxies, the 
catalogs appear to faithfully recover nearly all visible 
sources (see Figure [lj . We discuss this issue in more 
detail below. The field-of-view is a circle of 1.2° diame- 
ter. The pixel scale is 1.5 arcsec pixel -1 and the image 



point-spread function has a FWHM of 4.6 arcsec. Ex- 
posure times range from 541 to 4019 seconds per filter 
among the galaxies in our sample. The exposure times 
are equivalent in the two channels for all of our galax- 
ies. Quoted magnitudes are on the AB system and m = 
corresponds to 1.4 x 10~ 15 and 2.05 x 10~ 16 erg s _1 cm -2 
A^ 1 for the FUV and NUV bands, respectively. 

We use the NUV astrometry and aperture photometry 
to define the object sample and use the FUV matched 
aperture photometry to define the FUV — NUV color. 
Among the set of available apertures, we select the 7 pixel 
aperture as a compromise between including a significant 
fraction of the light from the source (7 pixels corresponds 
to approximately 2 FWHM) and minimizing contamina- 
tion from nearby sources. These apertures are plotted in 
Figure[T] We trim the full object catalog for each galaxy 
using a chosen outer image edge, beyond which the cata- 
logs appear to often have spurious sources, and an inside 
radius from the target galaxy, within which the catalogs 
are clearly deficient in detecting sources, perhaps due to 
source confusion. The outer cut is defined relative to the 
field-of-view center and occurs typically at a radius of 
about 1300 pixels. The inner cut is defined relative to 
the target galaxy nucleus and is typically at a radius of 
about 1.5 optical radii. The adopted values for each of 
these cuts are based on visual inspection of the images 
and source catalogs and are given in Table 1. An error 
in the placement of the inner cut will either cause a spu- 
rious drop in counts at small radii (if the cut is placed 
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Fig. 2. — Color-magnitude diagram of UV knots. The colors of 
two single s tellar population mo dels of age 100 and 360 Myr are 
shown from Bianchi ct al. (2005). 
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Fig. 3.— Correlation function of blue (FUV-NUV < 1) UV- 
selected knots for the combined final galaxy sample. 

at too small a radius), or preclude us from measuring 
the correlation function at small radii. An error in the 
outer radius can affect the knot counting statistics (if the 
choice does not properly exclude the area with spurious 
detections), but is clearly identifiable by a non-zero level 
of the background over a range of radii. 

Using the object catalogs, we calculate the angular cor- 
relation function, £(r), of sources. The correlation func- 
tion is defined in the standard manner — the excess prob- 
ability of finding a source at an angular distance r from 
the reference source. A random distribution of sources, 
such as that of the background field, will average out to 
a correlation function of zero at all radii (ignoring exotic 
effects such as gravitational lensing by the target galaxy 
or obscuration by galactic dust at large radii, for an ex- 
am ple of the effect s of the latter at optical wavelengths 
see IZaritskvl (|1994l )). 

Although a natural choice is to use all the objects in 



the catalog to construct £(r), judicious choices can en- 
hance the contrast of the extended disk UV knots to the 
background population. Previous studies in individual 
galaxies have stressed the blue UV c olor of the most evi- 
dent population of oute r disk sources (jThilker et al.ll2005l : 
iGil de Paz eiTall [2005) , so one advantageous possibility 
is to focus on the bluer of the detected sources. Figure 
[2] contains the color-magnitude distribution of sources 
in the NGC 628 field, which is one of the galaxies in 
the sample for which the extended disk UV knots are 
most evident. The Figure includes both blue and red 
plumes. The blue plume is where the previous extended 
disk sources were identified, and we have included la- 
bels indicating the corresponding color for single stel- 
la r populations at two ages drawn from the modeling 
of I Bianchi et al . (2005). When discussing results below, 
we will refer to the following populations: "all" sources 
have FUV - NUV < 10 and NUV < 25, "bright blue" 
sources have FUV - NUV < 1 and NUV < 23, "blue" 
sources have FUV - NUV < 1 and NUV < 25, "bright 
red" sources have FUV - NUV > 4 and NUV < 23, 
and "faint red" sources have FUV - NUV > 4 and 
23 < NUV. 

There are various choices that must be made when 
calculating £(r) from a sample of reference sources, as 
is the case here. First, the radial coordinate must be 
defined. The simplest choices are either purely angular 
(arcsec) or physical (kpc), but the galaxies are at dif- 
ferent distances and of different intrinsic sizes so neither 
is an optimal choice. Instead, we use a distance inde- 
pendent scaling radius: optical radii in angular units. 
With this choice, £(r) now describes the excess probabil- 
ity of finding a source at a certain number of optical radii 
from the central galaxy. Second, the density of source 
within different fields, and within a set number of square 
optical radii, must be considered. To account for field- 
to-ficld differences, we calculate £(r) in each field first, 
and then average the individual galaxy values of £(r) to 
obtain the final £(r). This approach has the apparent 
disadvantage that we are not using the entire sample to 
obtain a highly precise measurement of the mean back- 
ground, but instead measure the background from each 
individual image. However, it has the advantage that the 
measurement of £ (r) is not dominated by a field with an 
unusually high background. Within each image, we mea- 
sure the mean background using source counts at radii 
> 4 optical radii, which is the farthest that anyone has 
yet c laimed a detection of UV knots (jGil de Paz et alJ 
2005). Our initial size criteria have ensured that we have 
a sample of galaxies for which a single GALEX image has 
significant coverage beyond 4 optical radii from which we 
can measure a background source density. The final val- 
ues of £(r) are expressed as a surface density of objects 
in units of UV knots, or stellar cluster candidates, per 
square optical radius. 

3. RESULTS 

In Figure[3]we present £(r) for blue sources. There is a 
strong signal within two optical radii, and no significant 
excess of sources thereafter. The flat correlation function 
beyond two optical radii argues against certain poten- 
tial systematic problems and confirms that defining the 
mean background at radii > 4 optical radii is reasonable. 
The transition in £(r) is sharp and probably underrepre- 



4 



Zaritsky & Christlein 



10 



- \ 


i 




II 


- K / 

r \ / 






II 



1 2 3456789 

Radius [Optical Radii] 



Fig. 4. — Correlation function of UV-selected knots as a function 
of knot selection parameters: all sources (FUV — NUV < 10, 
NUV < 25) solid line; blue sample ( FUV - NUV < 1 and NUV 

< 25) dotted line; bright blue sample (FUV - NUV <1 and NUV 

< 23) short dashed line; bright red sample (FUV - NUV >4 and 
NUV < 23) long dashed line; and faint red sample (FUV - NUV 
> 1 and NUV> 23) dot-dash line. 
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Fig. 5. — Correlation function of UV-selected knots for the 
individual galaxies. Correlation functions from all sources in left 
panel, and only from blue (FUV-NUV < 1) in right panel. The 
correlation functions have been shifted vertically for clarity. The 
galaxies are in the same order as in Table [T] beginning with NGC 
628 at the top of each panel. 

sented here due to possible increasing incompleteness in 
the source catalogs toward smaller radii (high source den- 
sities). Some incompleteness is suggested by the negative 
£(r) values at small r for the two red samples of objects 
(Figure [4]), which might be expected to be background 
objects and therefore uniformly distributed. However, 
other explanations for the defic it at small r, s uch as ex- 
tinction of background sources (Zaritsky 1994), have not 
been excluded. 

As previously mentioned, the strength and clarity of 
the correlation function signal may depend on the color 
of the population. In Figure [4] we plot £(?') obtained 



using all objects and isolating blue, bright blue, bright 
red, and faint red sources. The strongest signal comes 
from using all the objects, and there is no signal from ei- 
ther red population. This combination of results suggests 
that there are some sources with intermediate colors that 
contribute signal. However, the stronger signal obtained 
inside of two optical radii using all sources is countered 
by a poorer baseline for £ (r) beyond two optical radii due 
to the poorer discrimination of background sources. For 
this reason, we focus on results using the blue sample. 

A key concern in any stacking analysis is whether the 
final result reflects a universal trend or is dominated by 
a small subsample of sources. In Figure [5] we plot the 
correlation functions for each of the 11 galaxies. One of 
the galaxies, NGC 1672, has no data at r < 2 optical 
radii and is excluded from further discussion here. Five 
galaxies show a significant rise in £(r) for small r (r < 2 
optical radii): NGC 628, NGC 4303, NGC 4736, M 83, 
and NGC 5474. We define what is significant here by cal- 
culating the number of knots detected in the "all" sample 
and requiring that this value be sufficiently large that it 
is statistically unlikely to happen in one galaxy in a sam- 
ple often. In fact, the least significant of these detections 
(NGC 4736) is likely to happen only 4% of the time by 
random chance (in other words, we have 96% confidence 
that our detections of knots in this galaxies is not due 
to random fluctuations). If we remove these five galaxies 
from the sample and combine the remainder, we find no 
significant positive correlation signal within two optical 
radii (Figure [6|) . The five galaxies for which there is no 
detected signal are farther on average than those that 
produce a signal (17.8 vs. 8.9 Mpc), suggesting that we 
may be less sensitive to the knots and more susceptible 
to source confusion, although excluding NGC 7479 from 
the calculation of the distances makes the difference less 
striking (13.5 vs. 8.9 Mpc). Using binomial statistics, we 
calculate the smallest intrinsic fraction of galaxies with 
knots for which there is at least a 10% chance of finding 5 
or more positive detections in a sample of 10. This value 
is 0.27, so we conclude that with 90% confidence at least 
27% of spirals have UV knots at radii between 1.25 and 

2 i?25- 

To examine the possibility that we are missing the UV 
knots in some galaxies due to their greater distances in 
more detail, we consider only blue sources (those with an 
age < 360 Myr). The typical uncertainty in the combined 
£ (r) for galaxies that show no evidence individually of ex- 
cess sources (Figure[6|) is ~ 1 in our units, which suggests 
that even if £ (r) was a factor of ~ 3 lower in these galax- 
ies we would have significantly detected a difference from 
zero in at least the innermost radial bin. The different 
mean distances of the samples corresponds to a factor of 
4 (2.3) in luminosity or 1.5 (0.9) mag, with or (without) 
NGC 7479. In Figured] we show that even if we raise the 
magnitude limit by 2 magnitudes (bright blue vs. blue 
samples) we retain a positive signal from the entire sam- 
ple. We conclude that distance alone is not responsible 
for the lack of a detection of extended disk UV-knots in 
the subsample of five galaxies. 

Finally, we examine whether we can enlarge the sample 
by relaxing some of the selection criteria. Our criteria, 
which are superposed on whatever criteria were applied 
by the original observers, are angular size, inclination, 
and morphological type. The quantitative cuts we adopt 
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for the first two are in some way arbitrary. For example, 
the need for face-on galaxies does not dictate one use an 
inclination value of 45°. Morphological type is the only 
one of the three for which a change affects the physical 
character of the sample, and therefore we do not relax 
that criteria. 

We begin by modestly modifying the inclination cut 
from 45° to 50°. This change increases our sample size 
from 11 to 17 (adding NGC 1961, 3351, 3368, 6902, 
7741, and 7793). Only one of these galaxies individu- 
ally show a significant rise in the correlation function at 
small radii (NGC 6902), where significant is defined as 
described previously, and the combined correlation func- 
tion decreases slightly so that the number of knots for the 
entire sample between 1.25 and 2 optical radii declines 
from 33±5.8 to 19.6±3.8. Although the numbers decline, 
perhaps because the mean distance of these galaxies is 
larger than that of the base sample, the significance of 
the detection of a source excess remains > 5<r. Increas- 
ing the inclination cut by an additional 10° only increases 
the sample by 50%. Large gains in sample size cannot be 
achieved by modest changes in the inclination cut, and 
the analysis becomes much more ambiguous due to the 
projection effects. 

We now modify the angular size criteria. The primary 
consideration here is that angular size maps onto distance 
and distance affects our analysis in two ways. First, the 
knots of the more distant galaxies are naturally fainter 
and hence the data do not reach as far down the lu- 
minosity function. This change reduces the number of 
knots detected significantly. Second, these fewer knots 
are embedded in a relatively more numerous background 
(because the intrinsically scarce bright knots are now 
embedded in the large background population of faint 
sources). Because selecting galaxies at larger distances 
reduces both the absolute numbers of detected knots and 
the contrast, we expect a sharp decrease in our ability to 
measure a positive correlation function. We have identi- 
fied a sample of 13 independent galaxies (i.e. not in ei- 
ther the base or inclination-extended sample) that have 
2 < a < 4 and are at distances < 50 Mpc. These were 
winnowed from a sample of 18 that satisfied the same 
criteria. The rejected galaxies either had bright or mul- 
tiple companions or there were large scale problems with 
the images (variable reddening or scattered light). At 50 
Mpc, the brightest of the knots observed in the nearer 
sample have an apparent magnitude of ~ 24.5. This sam- 
ple shows no significant excess (integrated counts from 
1.25 to 2 optical radii are 0.92 ± 1.65). To detect this 
value at 6cr, comparable to what is done with the base 
sample, would require a sample that is 115 times larger 
(i.e. 1500 galaxies). It is therefore not possible with the 
GALEX data set to exploit large numbers to offset the 
drop in detection sensitivity and relative contrast to the 
background. We use the results from the base sample in 
the discussion. 

4. DISCUSSION 

The identification of extended disks in at least a quar- 
ter of disk galaxies indicates that this is not a fringe 
phenomenon that affects some systems with rare specific 
histori es, such as recent interac ti ons (e .g . M 83 and NGC 
4625: iKobuinickv fc Skilmmnl (fl995h : iBush fc Wilcotsl 
(200i)). Instead, we conclude that spiral galaxy disks 
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Fig. 6. — Correlation function of UV-selected knots in a selected 
subsample of galaxies (circles with error bars). We removed the 
galaxies that show the strongest signature of extended disk sources 
(NGC 628, 4303, 4736, 5474, and M 83). The correlation function 
for the combined sample is the solid line and shown for comparison. 

are commonly at least 2 times larger than indicated by 
their optical radii and that in physical size they can ex- 
tend out to 26 kpc (2 x average optical radius, which 
is 12.9 kpc for this sample). While the average optical 
radius is 12.9 kpc, the galaxies in the sample radii span 
from 4.7 to 20.9 kpc. The five galaxies for which we de- 
tected correlated signal on an individual basis span the 
same range, suggesting that this phenomenon is not lim- 
ited to systems with abnormally large or small optical 
radii. 

Integrating the correlation function between radii of 
1.25 and 2 optical radii results in a value for the ex- 
pected excess number of sources of 75 ± 10 knots for 
galaxies with extended disks. This number must be 
treated with some caution because it applies to sources 
brighter than NUV = 25 for a set of galaxies over a 
range of distances. If between 5 and 10% of these clus- 
ters have associated Ha emission, as suggested in pre- 
liminary analysis of obse rvations of knots in NGC 628 
(|Herbert-Fort et al1l2007f) . then a long slit pointing of 
an edge-on galaxy should yield between three and seven 
knots. This expectation is in good agreement with the 
results of blind slit observati ons of edge on galaxies by 
IChristlein fc Zaritskvi (|2007al ). This ratio of Ha emit- 
ting sources to non-emitting sources is entirely reason- 
able given the age estimates of ~ 20 Myr for H II regions 
and ~ 400 for GALEX knots (see Figured]). Establish- 
ing the correspondence between the UV-knots and Ha 
sources is important because we will then be able to 
connect the kinematics obtained from the Ha observa- 
tions to the entire population of UV-knots. In particular, 
IChristlein fc Zaritskvi (|2007al) have established that the 
Ha sources follow the disk rotation curve out to 2i? 2 5 and 
hence lie in a dynamically stable, differentially rotating 
disk that continues from the inner disk. By inference, 
the UV knots are also then part of a dynamically cold, 
extended disk. 

Connecting the UV-knots to other components, for ex- 
ample H I or dust, will enable further study of the evo- 



6 



Zaritsky & Christlein 




2 3 456789 

Radius [Optical Radii] 




5 

FUV - NUV 



Fig. 7. — Correlation function of UV-selected knots in a selected 
subsample of galaxies that show the strongest signature of extended 
disk sources (NGC 628, 4303, 4736, 5474, and M 83). Solid line 
is the composite for the entire sample and open circles show the 
results excluding M 83. 

lution of material at these radii. iThilker et alj (2005) 
found a general correspondence between H I filamentary 
structure and UV- knots in M 83, suggesting that H I 
could be used to trace the extent of recent star forma- 
tion in extended disks. We find a correspondence be- 
tween the radial range over which we find UV knots 
and that ove r whic h lOO^tm emission was detected by 
iNelson et al.l (Il998f). In their most sensitive compos- 
ite sample. iNelson etall lfl998| ) found 100/im emission 
with 2cr confidence out to 2 optical radii, exactly the 
radius out to which we find UV knots. Furthermore, 
INelson et al.l (|1998f) show that the emission they de- 
tect comes from a disk- like component, again consistent 
with the geometry foun d here (in combination with t he 
kinematics presented bv lChristlein fc Zaritskvl (|2007aD ). 
Subsequent studies hav e found more direct evi dence for 
dust in extended disks (jPopescu fc Tuffs! I2003D and the 
correspondence of dust to these areas of star formation 
suggests that exotic transport mechanisms from the in- 
ner disk are unnecessary. 

To determine what the ages, masses, and contribution 
to the outer disk surface brig htness of these knots might 
be, we have run PEGASE.2 (jFioc fc Rocca-Volmeranjg 
119971 ) models of clusters. Rather than the standard sin- 
gle stellar population models, we vary the length of the 
burst (modeled by a Gaussian star formation rate in time 
of dispersion, S Myr). We include models with S = 1 
Myr and 100 Myr in Figure [51 These models are other- 
wise standard: Salpeter initial mass function (IMF) for 
0.1 < M < 12OM , and default PEGASE assumptions 
regarding winds, reprocessing of metals, nebular emis- 
sion, extinction. Changes in the initial metallicity or 
IMF do not produce changes that are sufficiently large 
to impact the conclusion discussed here. We show mod- 
els of clusters with total masses of 10 3 and 10 6 Mq in 
Figure[8j The mass of the 10 3 M Q cluster is only about 
a factor of two higher than the minimum mass cluster 
that is likely to have a star with M > 20M Q (the lowest 
mass star that will power an H II region), while the up- 



FlG. 8. — Models for the evolution of stellar clusters superposed 
on properties of observed UV knots in NGC 628 (corrected to a 
distance of 8.9 Mpc). Solid line represents a clusters with a mass 
of 1000 Mq with a Gaussian burst timescale of 100 Myr, while the 
dotted line represents the same mass cluster with a burst timescale 
of 1 Myr. The dot-dashed line represents a 10 6 Mq cluster with 
a burst timescale of 1 Myr. All models plotted for an assumed 
distance equal to the mean of the galaxies that show the excess 
knot population (8.9 Mpc). Models do not include any extinction. 

per mass limit appears to represent an upper bound to 
the likely cluster mass range given the dearth of observed 
clusters with NUV < 19. With this range of models, and 
some reddening variations, it appears entirely possible to 
have simulated clusters fill the entirety of the region in 
color-magnitude space occupied by the majority of de- 
tected sources. The two models with different S values 
demonstrate the sensitivity of the peak luminosity on the 
duration of the star formation episode, and signal the dif- 
ficulty in assigning unique ages or masses to candidates. 

We conclude from Figure [5] that the sources in the blue 
plume arc likely to be clusters younger than 400 Myr in 
the mass range 1 3 < M < 10 6 M Pl , in concurrence with 
earlier findings bvlThilker et alj (|2005l ); lGil de Paz et all 
(|2005h : iBianchi et all (|2005l ). If these clusters form uni- 
formly during those 400 Myr (consistent with the rough 
Ha statistics available) , and if 50% of spirals have these 
populations, then either 50% of spirals have such sources 
in their outer disks all the time, these sources are gen- 
erated on timescales of a few hundred Myr in spirals, or 
something in between. Either alternative, or a combina- 
tion, indicates that triggering star formation in the outer 
disk is common. If 50% of spirals have these all the time, 
and if they typically have 70 such knots between 1.25 and 
2 optical radii, then over the course of 10 Gyr they will 
have had 70(10/0.4) = 1750 knots. If these knots have 
a typical mass of 10 4 Mq, then the outer disk at these 
radii contains a stellar mass of ~ 1.7 x 10 7 Mq. Using 
the predicted evolution of the V-band magnitude from 
the models, we integrate (using broad, lGyr, time bins 
which results in a slight underestimate of the V magni- 
tude) to obtain the V magnitude of these stars currently. 
For an assumed distance of 8.9 Mpc we find V ~ 17.5. 
For a middle ground value of a of 4 arcmin, the area be- 
tween 1.25 and 2 optical radius equals 27596 sq. arcsec, 
implying an average surface brightness of ~ 29 mag/sq. 
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arcsec in V, which is comparable to what the deepest 
studies reach (although with the expected surface den- 
sity gradient expected from 1.25 to 2 optical radii, this 
population may be easily detectable at 1.25 optical radii 
and nearly impossible to detect near 2 optical radii) . We 
have ignored reddening due to the low inferred dust con- 
tent at these r adii from both extinction measurements 
(Zaritsky 1994) and far infrared luminosity (|Nelson et al.l 
1998). We conclude from our modeling that the observed 
rate of star formation in the outer disk integrated over 
a Hubble time does not grossly violate photometric ob- 
servations of disks, and that in some cases my be consis- 
tent with the obser ved optical surface brigh tness at large 
radii (jPohlen et al.ll2002t JPohlen k Truiilloll2006r ). These 
models are presented simply to begin a discussion of the 
characteristics of this population and clearly await both 
tighter observational constraints and more detailed and 
comprehensive modeling. 

5. CONCLUSIONS 

We have studied the distribution of UV-selected 
sources surrounding a sample of nearby, face-on spi- 
ral galaxies with the intent of quantifying some prop- 
erties of the previously discovered population of star 
forming regions well beyond the optical radii of 



these galaxies dFerguson et al.lll998; 


Thilker et al.ll2005; 


IGil de Paz et al. 2005: Bianchi et alj 


2005). We find the 



following: 

1) Using the cross-correlation function of selected sub- 
samples of UV knots and the primary galaxy, we detect 
an excess of sources in 50% of our sample for 1.25 < R < 
2 optical radii. From the integrated counts of knots in 
this radial range, the detection of knots at extended ra- 
dius as a general, although not ubiquitous, phenomenon 
of spiral galaxies is demonstrated at > 5cr significance. 
We conclude, with 90% confidence, that at least 27% of 
spirals have such UV source populations. 

2) In the combined data set (11 galaxies), we find no 
evidence for excess sources beyond 2 optical radii. Indi- 
vidually, only M 83, exhibits an e xcess of sources beyon d 
2 optical radii (in agreement with lThilker et a l. (2005)). 
Therefore, while such star formation is relatively com- 
mon interior to 2 optical radii, it is rare beyond. 

3) We find the majority of the correlation signal comes 
from blue sources (FUV — NUV < 1), but that there 
is detectable signal from redder sources (1 < FUV — 
NUV), indicating the presence either of some sources 
older than 400 Myr or reddened sources. If they are 
older than 400 Myr, then they need to be fairly massive 
(M > 10 5 M©). Because of the rapid fading of clusters at 
UV wavelengths, the prominence of the blue plume does 



not rule out older clusters, nor suggest the uniqueness of 
the current episode of cluster formation. 

4) The average number of excess blue sources between 
1.25 and 2 optical radii around spiral galaxies is 33.0±5.8 
(75 ± 10 for those galaxies that show a statistically sig- 
nificant excess individually). In combination with the 
numbers of observed Ha knots (a few in edge on obser- 
vations (IChristlein fc Zaritsky! l2007aft an d 10% of UV- 
selected sources (|Herbert-Fort et al.l l2007) in NGC 628), 
we conclude that the formation rate of these knots could 
be uniform over the last 0.5 Gyr. 

5) Integrating the number of such sources over the life- 
time of a galaxy, we obtain an estimate of the avearge 
V-band surface brightness of such sources of ~ 29 mag/ 
sq. arcsec. Given the expectation of a gradient in surface 
density, the surface brightness corresponding to these 
sources may be within relatively easy reach at the in- 
ner boundary of the range (1.25 optical radii) and nearly 
impossible at the outer edge (2 optical radii). The sim- 
ple calculation suggests that the e xtended disk s seen 
in diffuse light in man y galaxies (Poh len et al.l 120021 : 
Pohl en fc TruhiTol [200l could arise from a long-term 
population of evolved UV-knots. 

The outer, or extended, disk environment is well pop- 
ulated by recently (t < 400 Myr) formed stellar clusters 
in a large fraction of disk galaxies. As such, it opens 
several interesting avenues for study, ranging from the 
physics of star formation in metal-poor, low densit y en- 
vironments (see for example Boissi er" et al.l ()2006l )). to 
their use as kinematic tracers in the disk-halo interface. 
These clusters are a common feature of spiral disks and 
demonstrate that the stellar disks of galaxies are often 
twice as large as indicated by the central, high-surface 
brightness disk. 
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Disk Galaxy Sample 
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